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We recently found that Gamma-Ray Burst energies and luminosities, in their comoving frame, 
are remarkably similar This, coupled with the clustering of energetics once corrected for the 
collimation factor, suggests the possibility that all bursts, in their comoving frame, have the same 
peak energy (of the order of a few keV) and the same energetics of the prompt emission Ey (of 
the order of 2 x 10^^ erg). The large diversity of bursts energies is then due to the different bulk 
Lorentz factor Fq and jet aperture angle 0jet. We investigated, through a population synthesis code, 
what are the distributions of Fq and 0jet compatible with the observations. Both quantities must 
have preferred values, with log-normal best fitting distributions and (Fq) ^ 275 and (0jet) ~ 8.7°. 
Moreover, the peak values of the Fq and 0jet distributions must be related - 0jg;^Fo=const: the 
narrower the jet angle, the larger the bulk Lorentz factor We predict that ~6% of the bursts 
that point to us should not show any jet break in their afterglow light curve since they have 
sin 0jet < 1 /Fq. Finally, we estimate that the local rate of GRBs is ~0.3% of all local SNIb/c and 
^2.5% of local hypernovae, i.e. SNIb/c with broad absorption lines. 
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1. Introduction 

The spectral energy correlations in GRBs are still matter of hot debate. The isotropic equiv- 
alent energy ^iso of the prompt phase of long GRBs correlates with the rest frame peak Ep of the 
vFv spectrum [^: oc Ef^-^. A similar correlation (obeyed also by short events - [|lT|]) exists 



between the isotropic equivalent luminosity Liso and Ep [ ]24[ ] : Ep<x L 



rO.5 
^iso ■ 



If GRBs emit their radiation within a jet of opening angle ©jet, the true energy £'y~£'iso6ljgj can 
be estimated For ~30 GRBs with known Ojet, Ey is tightly correlated with Ep [^]. 

The presence of outliers of the Ep — Ei^o correlation [^], [pO|], | [23| ] [^] and the presence of pos- 
sible instrumental biases [Q], [|I9|], caution about the use of these coiTclations either for deepening 
into the physics of GRBs or for cosmological purposes. However, even if instrumental selection 
effects are present, it seems that they cannot produce the correlations we see [ p^ [jT5|]. More- 
over, a coiTclation between Ep and Liso is present within individual GRBs as a function of time 
0, [0] [0]. 

A new piece of information recently added to the puzzle is that the energetics in the comoving 
frame (i.e. E[^^, L[^^ and E'^^^,^ are similar for all GRBs [15]. For about 30 GRBs we [15| found 
that £'iso(Liso)°= Fq and LpOcFo; in the comoving frame £'j'j,jj~3.5x 10^^ erg, L[^^^5y.l(f^ erg s^' 



and £'p~6 keV (see [ 17 1 for a theoretical interpretation). These results suggest that the Ep —E\ 



p ^ISO 



and Ep — Liso congelations are a sequence of different Fq factors. 



The comoving true energy E'^ turns out to be ~ 2 x 10^^ erg. In [15] we argued that to have 



consistency between the Ep — Ey and the Ep — Ei^o coiTclations we need d^^J^o = constant. The dis- 
tribution of Fq is centered ai^ound Fo=65 (130) in the case of a wind (uniform) density distribution 
of the circum-burst medium. 

These new findings prompted us to explore the possibility that the Ep — Ey and the Ep — Ei^o 
correlations result from all bursts having the same comoving Ey and Ep but different Fq and Ojet- 
Specifically, we [0] ask whether dj^t and/or Fq have preferential values or not, and if there is a 
relation between them. To this aim we have performed extensive numerical simulations, along the 
guidelines explained below. 

2. Simulation set up 

Fig. I] shows the Ep — Ei^o plane. The black points are GRBs belonging to the complete Swift 
sample of [|22|]. The large black dot coiTcsponds to our main assumptions, i.e. all bursts, in the 
comoving frame, emit £''~ 2 x lO'^*^ erg at £'p~1.5 keV independent of their Fq. E' is smaller (2(7) 



than the mean value derived in [ 15], in order to be able to reproduce GRBs lying quite close to the 



1/3 

Ep oc E-J^ line. GRBs with different Fq would lie on the EpOcEy line, giving rise to the Ep — Ey 
relation. Then, by assuming a given aperture angle 0jet we can calculate £"150. The GRB will move 

to the right by the quantity [1/(1 — cos d^et)] if 6'jet>l/ro, and by the quantity 2Fq otherwise. In the 

1/3 

latter case, the relation between Ep and Ei^o becomes Ep <=< E^^^ . This implies that region (III) of 
Fig. |T] is forbidden. The other forbidden regions are region (II) because this would con^espond to 
0jet> 90°, and region (I) because we assume l<Fo<8000. All our simulated bursts will then lie 
on the white part of the plane. The distribution of the simulated bursts in this plane depends on the 
chosen distributions of Fq and ©jet- We thus have a tool to find what are the best fitting distributions. 
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Figure 1: Rest frame plane of ORB energetics. The large black dot corresponds to all bursts having the 
same E'^ and E'y. For a given To, the burst moves along the line Zspeak °= £y- Since we assume 1 < To < 8000, 
regions (I) are forbidden. Since all our simulated bursts have 0jet<9O°, they cannot lie in region (II). For 

1 /3 

small To, the beaming cone ~ 1 /Fq can become wider than 0jet. This introduces the limit fipeak °^ Ejo ^^'^ 
bursts cannot lie in region (III). Black dots correspond to the real GRBs of the Swift complete sample [p2[]. 



The steps are: i) select a redshift from the assumed redshift distribution (that is taken from 
[22], which includes an evolutionary term); ii) select a Fq and calculate Ep and Ey; iii) select a 6ljet 
and calculate ^iso; iv) chose a viewing angle and decide if it is pointing at us or not; v) calculate 
the peak flux in the appropriate band (assuming a typical Band spectrum) and decide if the burst 
belongs to the complete Swift sample [|^] or not. Bursts in this sample have a peak flux larger than 
2.6 ph cm^^ and almost 90% of them have a measured redshift. The steps are repeated until 
the number of simulated Swift bursts matches the real ones. Finally, we repeat 1,000 times each 
simulation to see how many times we can get a reasonable agreement with several observational 
constraints. First, we compare the simulated points of the complete Swift sample with the real 
ones in the Ep — Eiso plane. Then we compare them also in the observed planes £'p'"*-Fluence and 
£'p'"*-Peak Flux (irrespective if the redshift is known or not). Finally, we compare the distribution 
of simulated vs real flux and fluences of the BATSE and GBM bursts (down to limiting values that 
are not affected by incompleteness). 



2.1 Results 

We performed several simulations considering first that both Fq and 6ljet have no prefeiTcd 
values, i.e. assuming that they are distributed as power-laws, changing the coiTcsponding slopes. 
None of these cases is in agreement with the data. Then we assumed a broken power law either 
for Fo or for 0jet, or for both. For the latter case we do find some agreement, but the distribution 
of the simulated points in the Ep — ^iso plane describes a linear correlation, instead of the observed 
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Figure 2: Simulations assuming log-normal distributions of 0jet and To and the relation Qj^j ro=constant. 
Left panel: simulated points and real data (black) in the Ep — E\^o plane. Yellow points are all simulated 
bursts, blue points are those pointing at us, red contours are the distribution of simulated bursts (1 and a) 
brighter than the peak flux limit of the ^wi/r complete sample (i.e. 2.6 ph cm^^ s^'). Right panel: Simulated 
(contours) and real points (black squares) are compared in the £'p'"'-Fluence observational plane. 
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Table 1: Parameter values (/i and a) obtained by fitting a log-normal function to the distributions of Fq and 
0jet (Fig. 3), for all the simulated bursts (ALL), for those pointing to us (PO) and for those pointing to us 
and with a peak flux larger than 2.6 ph cm^^ s^' (the flux limit of the complete Swift sample) (PO Swift). 
For each distribution are re ported the three moments: the mode, the mean and the median. 



Ep oc E^^^. We then tried log-normal distributions both for Fq and 0jet. In addition we assumed that 
there is a relation between the average values of the two distributions. The best results are obtained 
with 6fjgj^ro=constant (Fig. ^). Note that the slope of the Ep — Ei^o correlation of bright bursts is 
harder than for faint ones (see the blue points in Fig. ^. But, curiously, these bright GRBs sample 
the distribution of the whole ensemble of bursts (yellow points) better than the fainter ones. This 
is because, if we improve our detector sensitivity, we preferentially see GRBs with larger opening 
angles. This makes them less energetic and enhances their probability to point at us. Fig. ^ shows 
(left panel) the distribution of Fq of all simulated bursts (black), those pointing at us (blue) and 
those (red) that are pointing at us and have a peak flux larger than 2.6 ph cm^^ ^-i g the flux 
limit of the complete Swift sample). The green points correspond to the few GRBs of measured Fq 
(left) or Qjet (right). Tab. |l] reports the parameters of the best fitting log-normal distributions values 
of Fo and d^ei for all bursts (ALL), for those pointing at us (PO) and for those pointing at us with 
peak flux larger than 2.6 ph cm^^ ^-i (t^g flux limit of the complete Swift sample) (PO Swift). 



4 



©jet and Tq ofGRBs 



G. Ghisellini^ 




10 



100 1000 



1000( 



10' 
10* 
lO' 



a 10 

•z, 

lO' 



10" 
10' 




10 



100 



[deg] 



Figure 3: Distribution of Fq (left) and of 0jet (right) of GRBs. Black circles: all simulated GRBs; open blue 
squares: all simulated GRBs pointing at us; open red circles: GRBs pointing at us with peak flux larger than 
2.6 ph cm^^ s^' (flux limit of the complete Swift sample); green triangles: the ^30 GRBs with Fq estimated 



from the onset of the afterglow 
on the right panel. 



on the left panel, and the 27 GRBs with measured Qjet collected in 



3. Conclusions 



The crucial assumption of this study is that all bursts have the same £'p=1.5 keV and 
2 X lO'*^ erg. Although there could be a dispersion of these values, our results still hold if the width 
of this dispersion is not larger than the dispersion of the observed quantities. The fact that these 
values are independent of Fq suggests that the dissipation mechanism giving rise to the prompt 
emission is not the transformation of bulk kinetic into random energy. If our assumption is true, 
then the — Ey relation is produced by the distribution of Fq values, and must be lineai^ (both 
Ep and Ey are proportional to Fq). In turn, the E^ — E{^o relation results from a distribution of jet 
aperture angles, with the caveat that, for small values of Fq, the radiation collimation angle is 1 /Fq, 
not ©jet- These bursts will never have a jet-break in the light curve of their afterglow, and could 
be mistaken as outliers. In our simulations we find that these should be about 6% of the GRBs 
pointing at us. Another important outcome of our study is that we can calculate the fraction of all 
GRBs (whether aligned or misaligned) with respect to SN Ibc, as a function of redshift. Taking the 



recent estimates of the SN Ibc of []18[], we find that, locally (i.e. up to z GRBs are 0.3% of all 
SN Ibc. 
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